Spatial reflection and renormalization group flow of quantum many-body systems 

with matrix product state representation 
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The property of quantum many-body systems under spatial reflection and tiie relevant physics 
of renormalization group (RG) procedure are revealed. By virtue of the matrix product state 
(MPS) representation, various attributes for translational invariant systems associated with spatial 
reflection are manifested. We demonstrate subsequently a conservation rule of the conjugative 
relation for reflectional MPS pairs under RG transformations and illustrate further the property of 
the fixed points of RG flows. Finally, we show that a similar rule exists with respect to the target 
states in the density matrix renormalization group algorithm. 
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Exploration of quantum many-body systems, particu- 
larly the translational invariant systems defined on lat- 
tices, is one of the most important topics in quantum 
physics and statistical physics. Nevertheless, to our ob- 
servation, the intrinsic attribute of the system under spa- 
tial reflection and the relevant physics have less been dis- 
closed so far. In particular, a question whether or not 
the species of quantum many-body states with matrix 
product construction possess inherently a reflection sym- 
metry [H is yet to be answered unambiguously. Here, 
we are motivated to reveal various categories of quan- 
tum many-body systems under spatial reflection and to 
explore the related property under the renormalization 
group (RG) procedure. 

The RG theory, including the seminal proposal of real- 
space renormalization by Wilson [2] and its renewed de- 
velopment of the density matrix renormalization group 
(DMRG) method [3], is one conceptual pillar of quan- 
tum many-body physics and particularly constitutes a 
key theoretical element to quantum critical phenomena. 
A theoretical picture of the standard DMRG algorithm 
could be formulated in terms of variational optimiza- 
tion within the representation of matrix product states 
(MPSs) UM- 

In fact, the generality of this mathemat- 
ical representation for quantum many-body states, in- 
corporating with the fact that the ground state of most 
quantum systems could be well approximated by a low- 
dimensional MPS, accounts unambiguously for the origin 
of the power of the DMRG algorithm. Recently, it was 
indicated that the MPS representation has a close con- 
nection with the concept developed in the field of quan- 
tum information, leading to significant progress, e.g., al- 
gorithms for periodic boundary conditions [5| , finite tem- 
perature 0, and simulating quantum systems of real 
time evolution Meanwhile, it was shown that a gen- 
eral RG procedure can be established upon the quan- 
tum state itself via MPS representation with properly 
defined coarse-graining transformations [§| . With respect 
to the Wilsonian RG scheme on Hamiltonians, this pro- 



posal suggests a specific rescaling approach to realize the 
scale separation for quantum many-body states. 

The main contribution of this paper are as follows. 
Firstly, by invoking the spatial reflection transformation, 
we show that apart from the symmetric states, the trans- 
lational invariant MPSs could have different attributes, 
that is, locally inequivalent to their reflectional counter- 
parts or differing from their reflectional counterparts only 
by local unitary transformations. Subsequently, we show 
that the conjugative relation of the reflectional MPS pair 
is preserved along the recurrent coarse-graining transfor- 
mations. Thus a rule on conservation of reflective relation 
for RG flows is indicated and the property of the leading 
fixed points is further investigated. Finally, we demon- 
strate elaborately that a similar law exists with respect 
to the target states in the numerical DMRG procedure 
regarding its variational nature of performance. 

Let us begin with the notation of the one-dimension 
translational invariant MPS: 



I*) = 



1 



d 

E 
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Tr(A^i...A^")|si,...,s^), (1) 



where the set oi DxD matrices {A", s = 1, • • • , d} param- 
eterize the iV-spin state with the dimension D < dJ^I"^ . 
The normalization factor is obtained as W = TrE^ , 
where E = X]s=i ® so-called transfer ma- 

trix with bar denoting complex conjugation. We now 
introduce a new state defined by a spatial refiection on 
I*), that is, \^rsi) = "^Jvl*) where is the parity 
operator for the iV-body system depicted by the action 
'Pn\s\ • • ■ Sat) = Isat • • • Si). In fact, for the present situa- 
tion with site- independent matrices {^^}, the reflectional 
counterpart state |^'r/i) is just an MPS represented by 
the set of matrices {A*^;} where A^.^^ = {A'^)'^ denotes 
the matrix transposition of A'^ . This can be easily seen 
from the equation 



I* 



rfl) 



1= y Tr{A' 



.A^")|sw,---,si) 



(2) 
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^—Y^Tr[{A^^f---{A^n^]\sir--,SN). 



Note that the transfer matrix of 
one of 1^) simply by 



rfl, 



is related to the 



Erfl 



Kfl = 



(3) 



which indicates that E and Erji have exactly the same 
spectrum. This leads clearly to the fact that any MPS 
has the same correlation length @,[n| with its reflectional 
counterpart. Furthermore, the overlap of two reflective 
MPSs is worked out to be 



7] = (*|* 



rfl) 



Tr{E^^) 



T2\N 



(4) 



TrE^ ' 

where E'^'^ = ® {^''V ■ Clearly, Eq. (g]) suggests a 

sufficient criterion for an MPS with reflection symmetry, 
that is, the specified matrix E'^'^ should have the same 
spectrum structure with that of the matrix E. 

The well-known multipartite states in quantum infor- 
mation, typically the Greenberger-Horne-Zeilinger state 
[iH, the cluster state [ll| and the MPS of Affleck- 
Kennedy-Lieb-Tasaki model [isj, are shown to be sym- 
metric under spatial reflection. Consider the cluster state 
as an example. By noting that the state has an MPS rep- 
resentation {A^ = {\ "), A^ = {]^ 0^)}^ one can work out 
that the matrix i?^^ has the same spectrum with that of 
i?, hence rj = \ according to Eq. More specifically, it 
is verified that E'^^ = {I(^X)E{I®X-'^), where the ma- 
trix X—[~'^ \) is an invertible transformation connecting 

the matrix and its transposition: {A'^)^ — XA^X^^. 

The translational invariant MPS without parity sym- 
metry, as will be shown below, exists in general. In fact, 
it is of interest to further distinguish two distinct cate- 
gories for the translational invariant states, namely, those 
locally inequivalent to their reflectional counterparts and 
those differing from their reflectional counterparts only 
by local unitary transformations. Specifically, let us con- 
sider a translational invariant MPS \'^) represented by 



A^ = 



1 

-I 



1 






1 9 



(5) 



It is shown that the corresponding transfer matrix E 
has eigenvalues — 1, — 1, (2 -|- 5^±\/4+^)/2 that are dis- 
tinctly different from those of the matrix i?^^ (except 
the case of g = 0). Therefore the specified reflec- 
tional MPS pair Ivf) and \'^rfi) are different according 
to Eq. (HI. Moreover, it can be shown that the two 
MPSs I^E*) and \^rfi) possess distinct correlation features 
[3], hence belong to different equivalence classes, i.e., 
I^'r/i) U (E) ■ ■ ■ ^U\'i/) where U stands for local unitary 
transformations. 

For the second example we consider an MPS j^*) with 



A' = 



1 
g 



1 




,A-' = 




9 



(6) 



Explicitly, the corresponding transfer matrix E has eigen- 
values 1 -I- g^, g, g, and the matrix E'^^ has 2g, 1, g^, 0, 
respectively. Therefore the MPS j^*) is different from 
its reflectional counterpart |^'r/i) in view that the over- 
lap between them is less than unity (apart from two ex- 
ception points of g = ±1). Interestingly, in this case 
the MPSs 1^) and \'^rfi) differ only by a local unitary 
transformation, i.e., \'^rfi) = U ^ ■ ■ ■ ^ f^l^) where 
U = |1)(1| + |2)(3| + |3)(2|. Definitely, the present ex- 
ample suggests a special category of translational invari- 
ant states that relate to their reflectional counterparts by 
non-trivial local unitary transformations. In general, the 
representative matrices of this special sort of MPSs sat- 
isfy (AO^ = UjiXA^X-^), where Uj is the represen- 
tative matrix accounting for the local unitary transfor- 
mation and X is an invertible matrix, say, it is obtained 
as X={^ °) for the present case of Eq. ([6]). 

Now let us consider the relation of the reflectional MPS 
pair under the RG transformation. Following Ref . , to 
perform the coarse- graining procedure for the state Iv?) 
in Eq. ([1]), one needs firstly merge the representative 
matrices for neighboring sites ^(p'^' = A^A'^. Then an 
appropriate representative for the equivalence class can 
be selected out via the singular value decomposition 



d' 



(7) 



where (pq) and (a/3) are understood as combined indices, 
and d' < imn{D^,cP} denotes the number of non-zero 
singular values of the matrix A'}^'^L . The state after one- 
step RG transformation can therefore be characterized 
by the new representative matrices 



AP ~^A- 



(8) 



Consider now the specified RG performance on the re- 
fiectional counterpart state \'^rfi) represented by {A^^^}. 
In view of the relation of the coarse-grained matrices 



A 



(pi) 

rfl 



^rfl^rfl 



(i(9p))T^ one has 



d' 



(^'■/0(a/3) 



1=1 



(9) 



Hence the RG transformation on \^rfi) gives rise to 



Kfi 



(10) 



Clearly, Eqs. ([S]) and pTI)l show that the relation of spa- 
tial reflection is preserved for the reflectional MPS pair 
under the RG transformation. In fact, the recurrent RG 
performance indicated by Eqs. ([7l)- (fT0|) suggests an in- 
triguing conjugative structure of RG flow for the trans- 
lational invariant states. This special flow conflguration 
will continue along the RG procedure until the states 



3 



reach their fixed points. Furthermore, since the corre- 
sponding transfer matrices after one-step RG transfor- 
mation are given by E' = E"^ and E'^^^ = [E^ f = {E'Y , 
the overlap of the reflectional MPS pair under the recur- 
rent RG performance is obtained exphcitly as 



TrE^ 



Tr[{E' 



TrE^ 



(11) 



where n = N/ x and we have denoted by E°° and E"^; 
{E°° Y' the transfer matrices of the two reflectional fixed 
points, i.e., E°° = lim^^oo E^. 

The above depicted conservation law of refiective re- 
lation for RG fiows is applicable for both the two cate- 
gory MPSs: those Y^) differing from \^rfi) by local uni- 
tary transformations and those \^) locally inequivalent to 
\^rfi)- For the former case, although and \^rfi) are 
viewed to be equivalent under the coarse-graining trans- 
formation, the attribute of spatial reflection is retained 
along the RG procedure even at the fixed point. In de- 
tail, let us examine the MPS of Eq. It is direct to 
calculate that the two reflectional fixed points are charac- 
terized by E^ = |$fl)($L| and E^^i = \^l){^rI where 
= (|00) + g^\ll))/{l + 9^) and \<^l) = |00) + |11). 
The corresponding representative matrices of fixed point 
MPSs |*°°) and are obtained respectively as 



1 






g 
















1 






g 



(12) 



and {{A'J^^Y = (A^)^, s 1, • • • ,4}. It is readily veri- 
fied that \^ffi) = U(g)- ■ ■(g)U\'i'°°) where the local unitary 
transformation U = + |4)(4| + |2)(3| + |3)(2|. 

For the situation specified by Eq. (O in which \'i>) 
and l^'r /i) possess different correlation feature, the cor- 
responding fixed points could be obtained similarly. In 
detail, since there is no degeneracy in the largest eigen- 
value of the transfer matrix E, the fixed points are char- 
acterized, up to an irrelevant normalization factor, by 
E°- = and E^fi = |$L)($fl|, where 

|1>k) = m + ^±^\U) {k^^4T7), (13) 



100) 



K + 1 — g- 



r(|01) + |10)) 



111). 



One can verify from Eq. ([4]) that for the corresponding 
fixed point states |(^'°°|*^,)| < 1. Notably, it turns out 
that the state differs from l^'J^j;) only by a local 

unitary transformation (15i] . Physically, it is understood 
that all correlation functions decay exponentially along 
the RG flow and become zero at the fixed point. There- 
fore the attribute of the fixed point described above is a 
general feature one exactly expects. 

So far, we have revealed various attributes and the rele- 
vant physics of RG fiows for translational invariant MPSs 
under spatial reflection. Now let us consider the DMRG 



scheme on the specified lattice system. Note that for a 
system without parity symmetry, i.e., Hrfi = 'P{H) ^ H, 
there is no reflection relation between the system block 
and the environment block in the DMRG algorithm any 
more. On the other hand, it is obvious that the systems 
H and Hrji have corresponding exact ground states re- 
lated by the reflection transformation. Hence it is in- 
teresting to explore whether the performance of DMRG 
algorithm could warrant the reflection relation between 
target states of reflective systems. 

In detail, let us look into the DMRG procedure with 
B • •B configuration for one-dimensional spin chains. 
The standard DMRG iterative performance could be de- 
scribed as follows. Suppose that the superblock com- 
prises two blocks and two spins in between at a certain 
step. The system block B^ contains spins 1, • • • M — 1, 
and the environment block B^ contains M+1, - ■ ■ , 2M — 
1. The states of two spins in between are denoted as \ sm) 
and Isjf), respectively. The target state, i.e., the ground 
state of the superblock has the following form 



I*) = 



E 



D 

E 

= 1 a, 13=1 



Ml 



(14) 

where the orthonormal bases \a) are obtained from 

previous steps and the tensor A'^^^"" is determined such 
that the target state minimizes the energy. From Eq. 
(fT4|) . the reduced density matrices of the left and right 
half superblock, B^* and •B^, are derived directly by 
virtue of the following singular value decomposition 



dxD 



' ,{s'^,..l3)i 



(15) 



where {sM,Oi) and {s'j^j,l3) are understood as combined 
indices and E is a diagonal matrix with elements (sin- 
gular values) Eq,/ , sorted in decreasing order, accounting 
for square roots of eigenvalues of the reduced density ma- 
trices. Then a truncation algorithm to achieve new sys- 
tem and environment blocks for the next step iteration, 
Bi» — !■ B^ and "B/j — s- B^, is performed by retaining 
only following D eigenvectors with the largest eigenvalues 



d D 



jj[M\,sm\ \L 

sjvf = l a=l 
d D 



M-1 



|sm), 



= E E^H"'"'"^^ 



A/-1 ^ PAf/ 



(16) 



SAf = l/3=l 



Here, C/]^^;^- and l/jfj^^- (a',/?' = 1, •••,/?) are 
just sub-unitary matrices truncated from ?7(sj,j,q),q' and 
^',(s'^j,/3)j respectively, and they fulfill the relation 

/ = ^([/[M],SM-)t[/[A/],SM 



4 



(17) 



In terms of the MPS representation, the target state in 
the above DMRG iteration procedure can be depicted 
distinctly as 

^YIM-1].sm+i . . . p.[l],«2M-l)|g^ . . .s2M-l), 



where 



we have 



used 



the notation |si • • • S2j\/_i) — 
|si • • • sm, s'j^ ■ ■ • S2M-1) and the summation indices {si} 
run over all the 2M spins. Note that the MPS here is site 
dependent, i.e., no longer translational invariant, and we 
have adopted periodic boundary conditions in Eq. (fT8|) . 
The DMRG procedure is now clearly phrased as that once 
the transformation A*"''*^ ([/[M].'**/ ^ I^M.s'm) jg ^jg, 

rived, then both system and environment blocks increase 
in length by one site and the algorithm is iterated until 
some desired final length is reached. 

The promising conjugated DMRG flow for the reflected 
system Hrfi is outlined below. It turns out that the 
target state of the superblock for the system Hr/i relates 
to the original one ^TE\i merely by an action of the parity 
operator |\E'r/i) =7^1^). Namely, one has 



1^ 



rfl) 



Tr(V"[^l''*^ • • • t/[*^~-'-l''**'^-iyl**'^'''Af (19) 



■ S2Af-l} 



where the tensors in the last expression are defined by 

Ylk],Sk _ ^Ylk],Sk-^T jj[k],S2M-k ^ (Jj[k],S2M-k'yr 



To demonstrate that the reflecting forms and 
of target states are preserved along the DMRG iteration 
for systems H and Hrfi, we need to prove (i) the formu- 
lated states (dH]) and ([TO]) minimize the energy of the two 
reflected systems simultaneously; (ii) the truncation al- 
gorithm of DMRG warrants that the resulted new repre- 
sentative matrices and target states satisfy repetitiously 
the indicated reflective relation. 

Point (i) is readily verifled since the expected values of 
Hamiltonians H and Hrfi over the states ^TE\\ and HI 
satisfy faithfully 



E = = {■i>rfl\Hrfl\^rfl)/W 



(21) 



with the normalization factor W = (*|*) = (*r/i|*r/i)- 
To demonstrate the point (ii), we note the following re- 
lation 



A 



a,f3 



{V EC/ )(sM,a)^{sl^,P) 
dxD 



/3',(sm,q)S/3'C/(s;j,/3),/3' 



(22) 



Consequently by virtue of the specified DMRG trunca- 
tion prescription one obtains the new representative ma- 



trices Vf 



l3',{sM,a) 



and U(. 



13', a — "aJJ' ^{slj,P),fS' ^ 

p 13' ~ ^ (}' p ' ^'^'^ corresponding recursive re- 
lations [cf. Eq. (dni] . This completes our proof that the 
presented forms and of target states are pre- 
served along the DMRG iterative procedure for the two 
reflective systems. 

In conclusion, we have disclosed the property of quan- 
tum many-body systems under spatial reflection and re- 
vealed a universal conjugative flow structure for both 
the RG scheme on translational invariant MPSs and the 
DMRG algorithm. An intriguing extension to high spa- 
tial dimensions via projected entangled pair states [l6j is 
awaited for us to explore further. 
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